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Abstract 

We study Schrodinger operators with periodic magnetic field in 
R^, in the case of irrational magnetic flux. Positive measure Cantor 
spectrum is generically expected in the presence of an electric poten- 
tial. We show that, even without electric potential, the spectrum has 
positive measure if the magnetic field is a perturbation of a constant 
one. 

Introduction 

Magnetic Schrodinger operators have been studied in Solid State Physics, 
especially in connection with the Quantum Hall effect, as well as on their 
own right. In a regular crystal 'physics' is periodic, i.e., the electric potential 
- caused by the background field of the ions - is a periodic function. Mag- 
netic fields - internal as well as external ones - are periodic as well, the latter 
ones typically being constant. Alas, as is well known, magnetic fields enter 
the Schrodinger operator through a vector potential, so that the resulting 
operator is not necessarily periodic. Indeed, it is so only in the simple and 
well-understood case of 'zero flux', where one has absolutely continuous spec- 
trum and band-structure ( [Birman fc Susliiia| , |1998| ; ^obolev| , |1999| ). Here, the 



magnetic flux (in units of flux quanta) is defined by 



* = 7^ / B{x,y)dxdy, (1) 
zn 



r 
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where B is the magnetic field and JF is a lattice cell (fundamental domain 
for the action of the group). Note that in our units with h = e = 1, the 
magnetic fiux quantum is just 27r. 

For integer or rational fiux the spectrum will still consist of bands (pos- 
sibly degenerating into points), but pure point spectrum is possible. 

For irrational fiux one expects Cantor spectrum (i.e. a nowhere dense set, 
no isolated points). The question is now: If B is constant, what 'defines 
the lattice', thus defining and the fiux? A potential V{x,y) = cos(27rx) + 
cos(27r?/) has periods (1, 1); one might consider any pair {m,n) of integers as 
periods of V, i.e., one can consider any coarser superlattice of Z x Z as the 
lattice of symmetry, but no finer lattice. Indeed, in this case one finds Cantor 
structure for a certain set of irrational values of $ ( [Helfi'er fc Sjostrand| , |1989D . 

On the other hand, V{x, y) = cos(27rx) has periods (1, c) for any real c; it 
does not define a fixed 'minimally coarse' lattice. Indeed, the Schrodinger op- 
erator with constant magnetic field and this potential V has band spectrum. 
This is still true for every potential V{x,y) = Vo{x) with reasonably non- 
degenerate Vq. Now, if we perturb such a V by a periodic Vi we expect Cantor 
spectrum, although this is not known. We only know that the Schodinger 
operator can be approximated in norm resolvent sense by operators from a 
natural algebra which have Cantor spectra ( pruber| , p001| ). 

Setting the Cantor issue aside, [Dinaburg et al.| ( |1997D showed that what 
survives under this kind of perturbation is the positive measure of the bands, 
although the bands might dissolve into a Cantor set. 

We ask the same question for Schrodinger operators with periodic mag- 
netic field B, without electric potential V. For constant B the spectrum is 
pure-point and infinitely degenerate (Landau levels). Is a periodic (zero-fiux) 
perturbation of B enough to cause the same effects as the potential V7 

In the course of giving an affirmative answer, we construct the generalized 
eigenfunctions and give estimates on the measure of the 'bands'. 



Outline 

In Section ^ we describe the setup and perform a first perturbation by one- 
dimensional magnetic potentials. For the general case, we give the direct inte- 
gral decomposition and express the operator in an appropriate basis 'moving 
along the fiber'. In Section ^ we give concise formulations of the main re- 
sults. The resulting double-infinite matrix problem is reduced in Section |^ 
to a problem which is almost diagonal in a sense made precise there. In 
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Section ^ we prove the main reduction lemma, using estimates on Weber- 
Hermite functions. Section |^ finishes the proof of the main results. 
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1 Setup and first perturbation 

For simplicity of the notation we assume that the lattice is Z x Z. Let B be 
an arbitrary smooth periodic magnetic field and $ its flux through a funda- 
mental cell of the lattice in units of flux quanta. Then we can decompose B 
as 

B = Bc + B^ with 
Be = 2n^ and 
Bz = B — Be 

Note that Bz has zero flux! Therefore, Bz has a periodic vector potential Az, 
and we can choose it to be of the form 

with periodic smooth and A^, where we introduced parameters eQ,ei for 
later convenience. 

Be is obviously constant, and we choose a vector potential 

My) = Be (i) (3) 



^0^ 

for it. In this gauge, the magnetic Schrodinger operator takes on the form 
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- Bey - .oA°(,))' + (i^ - e.A^x^y)^ 



(4) 



Note that, in order to emphasize the dependencies, we write out the argu- 
ments of A^ and A^ even though they are multiplication operators; in the 
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following we will keep sacrificing strictness of notation for better clarity in 
cases such as this one. 

First we will set ei = and study the effect of switching on Eq. The 
operator (|^) does not depend on x if ei = 0. Using Fourier transform on 
L^(Ma;), we can decompose L^{M?) into a direct integral L'^{Wy)d^ such 
that 

L\R') 3f^f, f^{y) = [ f{x,y)e-'-'^^dx, (5a) 
f{x,y)= [ /e(y)e2-«-rfe; (5b) 



H= H^di, H^ = V^{y)-\f^ with (5c) 

dy 

V^{y) = \{27x^-B,y-e,A\y)f . (5d) 

In the case Sq = £i = Q we are dealing with the Landau Hamiltonian. We 
will go through its analysis since we will use its eigenf unctions later on as a 
basis. 



1.1 So = = El 

V^{y) = \ {2ni - B^yf = f {y - (6) 

is a harmonic oscillator potential with frequency Be, shifted by (3^ with (3 = 
^ = 1 (we assume Be ^ 0, or else there is not much to do). If we denote by 

the Weber-Hermite functions, i.e. the standard normalized eigenfunctions of 
the harmonic oscillator with frequency 1, then 

^C,m= V^fim(/B;(2/-/30), mGZ+, (8) 

are the eigenfunctions of H^, and the corresponding eigenvalues are i?c (| + ^) 
Since the spectrum of is independent of ^ it coincides with the spectrum of 
if as a set, and both are pure-point. On the other hand, since H is invariant 
in X, it has infinitely degenerate spectrum. 
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1.2 £0 ^ = £1 



Introducing the shifted variable y = y — we can write the potential 

as 

B] 



niy) = ^ {y + eoA'' iy + (30y 



(9) 



Since is periodic and smooth (and therefore bounded), this is just a per- 
turbation of the harmonic oscillator potential. V^{y) will still tend to infin- 
ity as \y\ does, so that has discrete spectrum. Some simple estimates 
using test functions show that for small Eq the eigenvalues will be within 
Cm^o Hiax l^^l of the Landau levels, for some constants Cm (involving the 
maximum of \yflrn,{y)\) depending on m and Be only; in particular, they are 
independent of ^. A closer investigation shows that Cm is indeed bounded 
with respect to m. 

Note that is periodic in ^ with period §^ = ^- Therefore, the spectrum 

^ I 1 

of H = jjg consists of bands whose size is bounded by 2Cm£o max 1^4 |. 

These bands might degenerate into points. 



1.3 £0 7^ 7^ £1 

Now H is not independent of x any more. But at least it will be periodic 
in X with period 1, since A^ is periodic. Using Fourier series on L'^{[0, 1]^), 
we can decompose L^(]R^) into a direct integral -^^([0, l]x x My) such 
that 



H 



[0,1] 



f{x,y) 



[0,1] 



(lOa) 
(10b) 



1^ 

^ dx 



+ 27^^-B,y-eoA'>{y)) + 



+ I i^-eiA\x,y) 



;ioc) 



acting on functions periodic in x. Note that we keep denoting the fibre 
operator for the new direct integral, in order to avoid an inflation of 
notation. 
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If we now choose the basis (e^™^)^^^ in L'^{[0, l]x), defining an isomor- 
phism with £^(Z), and combine this with the isomorphism defining the direct 
integral above, we arrive at the direct integral L'^{M?) = £^(Z)(8)L^(]R) rf^ 
with 



L\R')3f^f, fUy)= f V/(x + /,i/)e-2-[«'+«+")^lc/a;, (11a) 
fix, y)= [ J2 4n(l/)e^^^(«+")^ d^; (lib) 
H = [ H^di, H^ = ]-\{2Ti{u + i)-B,y-e^A\y)Y + 

"'[0,11 ^ ^ 



[no 



Here, u is the operator of multiplication on £^(Z), i.e. {vg){n) = ng{n), and 
A^{y)ic is convolution with the Fourier series of A-^ix, y) in x: 

{A^{y)^g){n) = J^dilWn.iiy), (12a) 
A\,iy)= [ A\x,y)e~''-'^^dx (12b) 

J [0,1] 

Note that, of course, our basis functions belong to the domain of the operator. 

As a final step, we choose a special basis in -L^(M), namely the eigenfunc- 
tions (^c,m)mez+ described in subsection |L1], equation (|). Thus we arrive 
at 

Lemma 1. There is a decomposition o/L^(]R^) into Jj®^^ £^(Z x Z+) rf^ with 

L'{R') 3f^f, (13a) 
4n,-= / 5^/(x + /,y)e-2-^[«'+(«+")^lvl>^+„,^(l/)tixrfi/, (13b) 

^[0,1] xR ^^2; 

f{x,y)= [ A,n,^e2-^(«+")^vl/^+„,^(i/)rfe (13c) 
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such that H decomposes as 



J [0,1] 



Ln 



1 n V — -\ ^ — \ (n—k) 



Lk 



where 



l€Z+ fcGZ 

iez+ kez 



(14a) 
(14b) 
(14c) 

(14d) 
(14e) 
(14f) 



AlJp) = Bi j A\y + p)yni (^,y) (^/B^y) dy, 



(15a) 
(15b) 



— (fc) 



— (k) 



% jjj^\{y + p)ni (yW^iy + 13k)) (/B^i/) dy, 



\b, ^h{y + P) [^'i (/b; {y + 13k)) a„ (/b;^ 
- ni(^,{y + Pk))n'^(^/B',y)\ dy. 



(15c) 



(15d) 



As above, A^k denotes the k-th Fourier coefficient of A^ with respect to x, 
and analogously for (A^)"^. 

Proof. We have 

( ) = zJ zJ ^i,rn,n;l,kfu,k 

1&+ kez 

so that we just have to compute the matrix elements in the given basis, for 



all the terms in (11c). 



( |14b| ) These are just the Landau eigenvalues in the case Eq = = ei. 



(T^) & (15a) For the terms with coefficient Eq we have to compute the 
matrix element of the square term (A^)^ which is 

^/K Jj,Ay{y)n^ (/b: (y - + n))) (/b: - + n))) dy 



(14d) &; ( |15b| ) The terms with coefficient in QllcD give the matrix element 
of the mixed term {Bey — 27r(z/ + C,))A^, the calculation is the same as above. 



(|14e| ) &Z, ( |15c|) The coefficient £^ singles out the square term -k A^-k = 

[A^y^k, i.e. convolution with the Fourier series of (A^)'^. Its matrix element 
is 

/b: ^ (A^Vfc(2/)^- (/b: - + n))) (/b: - /^(e + fc))) dy 

so that shifting ?/ as above gives the desired result. 

dm & (|15|) The term with coefficient Ei is z o * ^ . The 
first part can be written -^oA^-k = A^ -k+A^-k-^ which has matrix element 



'Br 



A\^M^n. (^/B,{y-m + n)))nA^/B,{y-m + k)) 



+ A\.k{y)^l^ (j/ _ + n)) fi; /b; (2/ - + A:)) 



dy. 



On the other hand, we can use partial integration for the matrix element of 
the second part, which is A^ -k 



Br 



A^ 



'Br 



n^k{ynm [VBc {y - m + n))) n\ i^^/Br {y - m + k))) dy 

-A\_,{y)nrr. (^/B'riy-Pi^ + n))) nJ^/W^iy- + k)) 



B, A\.kiy)n'^ /b: (y - + n)) U^J /b: (y - + fc)) dy. 



These two parts add up to the desired result. 



□ 

~ ~(fc) (fc) 

Remark 1. Note that the functions A^i^m{p), ^imip)^ ^^imip)^ ^^imip) have 

period 1 in p. Also, due to the decay of the Weber- Hermite functions flm{y) 

~(fc) — (k) 

in y, the functions A^i^{p), A^i^{p) exhibit exponential decay in k. 



Remark 2. If depends on x only then ^^^^(p) = SkfiSi^mi^^y with the 
average {A^y of (A^)"^ with respect to x. If we assume that A^ = (we can 

always add a constant to A^ to achieve this, without changing the magnetic 

(f^\ 

field) then A\ „^{p) = 0. 

Remark 3. In the simplest non-trivial case A^{y) = cos{27iy), using parity of 
the Weber-Hermite functions we get 

where ^ = if m -|- Hs even, and 6;^^ = if m -|- / is odd. Similarly, 

^°i,m{p) = |5«,m + Q,„cos(47rp) + ci/,^ sin(47rp) , 

where ci^m = if m + / is odd, and di^m = if m + / is even. 

Remark 4. Creation and annihilation operators on harmonic oscillator func- 
tions yield the relations 



(/j) 

Using these we can express A'^i^m and A'^^^ solely in terms of A^, A^ and 
Weber-Hermite functions, without referring to their derivatives or multipli- 
cation by y. 



2 Main results 

Let P°„^ denote the projection on the eigenspace of the m-th Landau level. 
Lemma |l| tells us that the action of in this eigenspace, i.e. the part of 
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which is 'diagonal in m', amounts to: 

= (b, (i + m) + IclAl^^^ {I3{(, + n)) - (/3(5 + n))) /4.„,„ (17) 



hlA^'^^ m + n))- ei^ 1""'^ m + n)) 



,m,k 



For fixed ^ and m this is a one-dimensional difference operator with quasiperi- 
odic coefficients and exponentially decaying off-diagonal (i.e. k ^ 0) terms. If 
we choose to be independent of y as in Remark ^ there are no off-diagonal 
terms at all. If we furthermore choose A^{y) = cos{27iy) as in Remark]^ then 
([T7|) will be similar to the Almost Matthieu operator, with a slightly more 
comphcated potential. Indeed, if we look at terms of order up to Eq only 
it will be exactly the Almost Matthieu operator. In the case of constant 
magnetic field this observation goes back to [Hofstadt"er| ( |1976| ). 



If £0 and El are small enough then H will have invariant subspaces Em 
which are close to the eigenspaces of the Landau levels. We will construct a 
unitary transformation to achieve the following: 

Theorem 1. For every small enough Eo and Ei there is an M{eo^ei) such 
that M{eo^ei) ^ oo as both eo,Ei oo, and such that for m < M{eq,ei) 
the invariant subspace Em of H and the restriction of H to Em have a direct 
integral decompostion 

Em= r fmd^, (18a) 

J [0,1] 



Hri 



■.= H\e^= / Hm,^dC (18b) 

J [0,1] 

Furthermore, i?m,^ octs on g & £^(Z) as a one- dimensional difference operator 
with exponentially decaying coefficients close to those of {^n\): 

{Hm,^g) (n) = dm (/9(e + n)) g{n) + ^ a^(n - /c, p{i + n))g{k), (19a) 

fcez 

\\dm{-) - [Be {\ + m)- £o^^,„(-)] 11^2(51) < Co^o. (19b) 

^ \\am{k,-)e^\''\\\^^^^^^<CiEQ (19c) 

A:6Z\{0} 

for some 5 > 
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We denote by || ■ ||c2(si) the sum of the supremum norms of the derivatives 
of order 0, 1, 2. 

As in ( pinaburg et al.| , |1997| ) we can now make use of the resuhs of 



Dinabur^ ( |1997|) on ergodic famihes of operators. In order to apply these 



resuhs we need the foUowing assumptions: 

diophantine There are C > 0, k > such that > C/\n\'^ for all 

n G Z \ 0. Here, /? = = as before, and {■} denotes the fractional 
part. 

smoothness and y) are smooth: furthermore, all derivatives 
are analytic in < 5 for some common 5 > 0. 

Morse m is a Morse function on with exactly two critical points. 

Theorem 2. Let j3 be diophantine and A^, A^ smooth as defined above, and 
M > such that m ^-^ « Morse function with exactly two critical points 
for all m < M. Then there are eo,ei > depending on M such that for all 
Eq < eQ,ei < 6i the following are true for all m < M: 

1. There are 1-periodic measurable functions \m such that for every n E'L 
and almost every G [0, 1], Xm{/3{(, + ^)) is an eigenvalue of Hm,s^ and 
therefore in the spectrum of Hm- Furthermore, 

||A™(-)- [5e(| + m) -£o^^,^(-)]Loo(5i) <const£2. (20) 

2. There are 1-periodic measurable functions fm,i,kJ £ 1j+,n G Z, decay- 
ing exponentially in n such that for almost every ^ G [0, 1], 

ElW/^OI(^' + l)e'""'<oo, (21) 



and for every G Z, the series 

$„,5,fc(x, y)=Y. fn,i,n-km + A;))e2-^(«+")^'v^5+„Xl/) (22) 



ie2 



and all its derivatives converge uniformly in x,y. $m,^,fc is an eigen- 
function of Hm,^ and therefore a generalized eigenfunction of H, with 
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eigenvalue Xm{(3{C, + k)). Moreover, for every N > and for eQ,ei 
small enough (depending on N), 

const , , 

|$„^,C,;.(x,i/)|<^5^ (23) 

with the constant depending on ^, k, m, N. 

3. Hm is uniformly e^-close to hand structure: is unitarily equivalent 
to multiplication by the function Am(/5-). The Lebesgue measure of the 



spectrum of Hm is eo \ ranA^^| + 0{e, 



Remark 5. One can reduce the smoothness requirements somewhat (to an- 
alyticity in x of a finite number of derivatives in y), thereby weakening the 
resuh on decay of the generahzed eigenfunctions in y. 

Remark 6. Of course one can also include an electric potential into the pic- 
ture; Pinaburg et al.| ( |1997| ) did so in the case of constant magnetic field, i.e. 



A'^ = = 0. The point in our work is that the magnetic field perturbation 
alone is strong enough to deform the Landau levels into a spectral set with 
positive measure. 



3 Reduction 

is a double matrix operator on £^(Z_|_ x Z) with indices {m,n) G Z+ x Z. 
We decompose it as = Vi + A4i + Oi, where Vi is diagonal in m and n, 
0\ is off-diagonal in m and contains only the first row and the first column, 
and A^i is the remainder. Note that both 0\ and M.^ are of order Eq (we 
always assume e\ < Eq). 

Our goal is to find a unitary transformation U which kills the terms in O 
(they represent the interaction between different Landau levels). U should 
leave the rest of the structure basically untouched. We will show how to 
accomplish this for the interaction between the 0-th and all other bands; 
extending this to M off-diagonal terms is a trivial generalisation. 

The strategy is as follows: We define U as U = HjeN "where Uj elimi- 
nates off-diagonal terms up to (and including) order ^q. Each transformation 
is of the form f/,- = e*'^-' for a Hermitian bounded whose coefficients are 
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of order e^. We use the Baker-Hausdorff formula 



[B,A]o = B, 
[B,AU^ = [[B,Al,A] 



(24a) 

(24b) 
(24c) 



in order to find Wj. In fact, in the j-th step we will only have to consider the 
terms up to order Eq in this formula, which are the terms r = and r = 1. 
This gives us the equation 

O, + ^[V,,W^]=0, 
which reads as follows for the coefficients: 



Wi{m,n]l,k) = for m > 0, 
Oi{0,n]l,k) +t[Di{0,n;0,n) -Vi{l,k;l,k)]Wi{0,n;l,k) = 



(25a) 
(25b) 



Because of 



\Vi{0,n;0,n) -Vi{l,k;l,k)\ > 



n 1 — ^c"^ 



2^1 



A 



(0) 



AO 



1,1 



0,0 



A 



(0) 
LI 



— ^0 (||^°0,0 

-(0) 



\A' 



Ll\ 



El 



A' 



0,0 



A 



-(0) 
Ll 



> \B,l 



we can choose Eq small enough so that for all / > there is no small denom- 
inator problem when solving equation ( p5b| ) for Wi. 

Defining H2 = UlHUi and repeating the above steps we arrive at the 
following: 



Lemma 2. We define Wj inductively by [T>j,Wj] = lOj, and furth 
Uj = e*^^ Hj+i = U*HjUj. Then we have: 

closeness to diagonal operator 

\\Vj{m,n;m,n) - [B^d + m] - ^o^^,™] || 0-2(51) < ^o^j 



ermore 



(26) 
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off-diagonal smallness 



5^||O,(0,n;/,n)||^,(5,)(r + l)<7, (27a) 
J2 E mO,n;l,m^,^,,^e''\-'\r + l)<e,^, (27b) 



smallness of mixed terms 



E \\M,{m,n-l,n)\y^s^^{l' + l) < (m^+^ + 1)5,- (28a) 



E E l|A<.(rn,n;/,fc)||^.(^i)e2^l'^-'^l(f + 1) < (m^^^ + l)£i5, (28b) 
kez\{n} iez+ 



Furthermore, we have the relations 

7j+i = const 7j5j, (29a) 
Sj+i = Sj{l + const 7j), (29b) 

W,-(0, n; k) ~ ^^'^"'""'^'^^ as Z ^ oo, (29c) 

where the constants depend on s. 

Later we will see that 

6j < consts^o 
jj < (consts EoY 
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4 Estimates 

In order to prove Lemma ^ we only need to check the induction hypothesis. 
For j = 1 we have the following non-vanishing matrix elements: 

Vi{m,n;m,n) = + m) + ^elAS>m,m{P{C + n)) - eoAl^^(p(C + n)) 

+ lelA^t]m m + n)) - m + ^)) ; (30a) 

MM, n; I, k) \elA^^Z'^ m + n)) - e^A^^^'^ m + n)) , (30b) 
Mi{m,n-l,n) = \elA^^,i{(3{i + n)) - eoAl^^i{l3{i + n)) 

+ IslA^'-Z (/^(^ + ^)) - ^1^^!!!^ (/^(^ + ^)) ; (30c) 
Oi(m,n;0,n) = \elA^^^^{i5{i^ n)) - e^A'^^^^^ii^i^ n)) 

+ \e\A^'t,^ + n)) - !;|o m + n)) , (30d) 

Oi(m, n; 0, k) ie^A^I^,;'^ + n)) - eJ^'^Z^^ m + n)) (30e) 

— ^(0) 

Note that we may assume -A^^n ; = 0; we can always achieve this by changing 

~{0) 

A by a function of y only, this does not change B. A^^ ^ = cannot be made 

disappear in the same way because it involves the average of (A^)'^ instead 
of A\ 

If we rewrite the induction hypotheses (p6D-( |29^ ) in terms of the coeffi- 
cients of the vector potential we get: 

Yl \\l^oA^m,i m + n)) - Al, + n)) 

+ -^CiA'Z (/?(« + n)) - ^ + n)) < const (m'+l + l) (31a) 



EE 



< const (' 



1) (31b) 



Since the vector potentials A^, A^ are periodic in y, equations ( |31a| ) and ( pib| ) 
can be derived from the smoothness (analyticity) assumption and estimates 
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for 

Jr 

in 7]. For this we can use the estimates from the case of constant magnetic 
field ( pinaburg et al.| , |1997| , lemmata 6-8). Because of our Remark ^ there 



are no new estimates to prove. 



5 Finishing the Proof 

The columns of the tranformation U define new basis vectors b^^n,m, n ^ 
Z, m G Z+. From the induction hypothesis ( |27b| ) we get 

\k,n,m\{m^^^ + I)e2^l"-Jl < const . (32) 



This implies that for each m, we can apply the main results of [Pinaburg 



(|1997| ) to the operator H^, which is an ergodic family (in ^) of difference 



operators with exponentially decaying almost periodic coefficients. As a con- 
sequence we get part |I] and the eigenfunction expansion in part 0, equations 
dll) and (H), of Theorem |. 

Uniform convergence of ( |2^ ) and its derivatives follows from ( |52| ) and the 
relation (|16a|) for the derivatives. As in ( pinaburg et al.| , p.997| ), (^) also 



gives the estimate on the rate of decay in y, 

Finally, from part ^ we know that, for almost every C,, the spectrum of 
Hm,^ is given by {Xm{P{C, + ^)))nez- E^^rn,n are the corresponding eigenspa- 
ces, then Em,n '■= J E^,ni,nd^ are invariant subspaces such that the restric- 
tion of H to Em,n is unitarily equivalent to multiplication with the function 
^ XmiPiC, + n)). On the other hand, Hm = ®n&z^rn,n- Therefore, since 
Hm = J Hm,i d^i the operator Hm is unitarily equivalent to multiplication 
with the function Am- 

This together with the quantitive estimate ( pOD finally gives the assertion 
in part ^ about the measure of the spectrum. 
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